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Table 1. Uncertainties of parameters, Bulge (I = 15, ip = 225°, 5 hours) 
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Table 1 — Continued 
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Table 2. Uncertainties of parameters, LMC (J — 15, arbitrary ip, 5 hours) 
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££ ■ ABSTRACT 

Astrometric measurements of microlensing events can in principle 
P-h ! determine both the "parallax" fg and the "proper motion" \x of an 

OO \ individual event which (combined with the Einstein time scale ig) in 

,_i ! turn yield the mass, distance, and transverse velocity of the lens. We 

show, however, that the parallax measurements are generically several 
t— i \ orders of magnitude less precise than the proper-motion measurements. 

Fortunately, astrometric measurements by the Space Interferometry 
Mission (SIM) are simultaneously •photometric measurements, and 
q^ \ since SIM will be in solar orbit, these allow SIM to be used as a 

"^ | classical (photometric) parallax satellite. We show that SIM photometric 

£p-( parallaxes are of comparable precision to its astrometric proper-motion 

measurements. For / = 15 bulge stars, complete solutions with ~ 5% 
^ | accuracy in mass, distance, and transverse velocity can be obtained from 

.. ', about 5 hours of observation, 100 to 10,000 times shorter than would be 

. £h ! required for a purely astrometric solution of similar precision. Thus, it 

rS \ should be possible to directly measure the mass functions of both the 

bulge and the inner disk (including both dark and luminous objects) 
with only a few hundred hours of SIM observations. 

Subject headings: astrometry - Galaxy: stellar content - 
gravitational lensing - Magellanic Clouds 



1. Introduction 

Boden, Shao & Van Buren (1998) have shown that it is in principle possible 
to obtain complete solutions for microlensing events from a series of astrometric 
measurements using the Space Interferometry Mission or possibly ground-based 
interferometers. This would be extremely important if practical because two major 
questions that are difficult to answer on the basis of present-day data could then be 
easily resolved. 

First, after almost a decade of observations, the nature of the events currently 
being detected toward the Large Magellanic Cloud (LMC) by the MACHO (Alcock 
et al. 1997b) and EROS (Aubourg et al. 1993) collaborations are a complete mystery. 
On the one hand, the observed optical depth r~2x 10~ 7 is an order of magnitude 
higher than expected from known populations of stars. On the other hand, if the 
lenses lie in the Galactic halo and so comprise of order half the dark matter, then 
their masses (inferred from the event time scales and kinematic models of the halo) 
are of order half a solar mass. Thus, the objects could not be made of hydrogen 
or they would have easily been discovered from star counts (Alcock et al. 1997b 
and references therein). Direct measurements of the mass and distance of the lenses 
would unambigously resolve this question. 

Second, there appears to be a large excess of short time scale events toward the 
Galactic bulge relative to what would be expected if bulge stars had a mass function 
similar to that seen in the solar neighborhood (Han & Gould 1996). The events 
would easily be explained if the mass function were more steeply rising toward low 
masses (Zhao, Spergel, & Rich 1995; Han 1997), but recent observations of the bulge 
by Holtzman et al. (1998) show that the bulge luminosity function is very similar to 
the local one. By now hundreds of events have been discovered toward the bulge 
although only about 50 have been published (Udalski et al. 1994; Alcock et al. 
1997a). If individual masses, positions, and velocities of even 10% of these could be 
measured, our knowledge of the bulge population (both dark and luminous) would 
be dramatically increased. 

In addition, the PLANET (Albrow et al. 1999) and MPS (Rhie et al. 1999) 
collaborations are currently searching for planetary systems by closely monitoring 
ongoing microlensing events seen toward the Galactic bulge. Ordinarily, these 
observations can yield only the planet/star mass ratio and their projected separation 
in units of the Einstein radius of the lens (Mao & Paczyhski 1991; Gould & Loeb 



1992). Complete solutions of the event would enable one to translate these quantities 
into planet masses and physical projected separations. 

At present, the only quantity routinely measured for all events is the Einstein 
time scale, t E , which is a complicated combination of the physical parameters that 
one would like to know, 

t E = — , ^ = TT I 1 ) 

where v is the transverse speed of the lens relative to the observer-source line of 
sight, fi is the proper motion, and 6 E is the angular Einstein radius, 

_(4GM\W D ol D os 

Here, M is the mass of the lens, and D \, D os , and D\ s are the distances between 
the observer, lens, and source. There are numerous ideas on how to get additional 
information about individual events, but these often require special circumstances. 
For example, if the source crosses a caustic in the lens geometry, then it is possible to 
measure the proper motion /i, and so from equation ([[]), the Einstein radius (Gould 
1994a; Nemiroff & Wickramasinghe 1994; Witt & Mao 1994). In fact, a variant 
of this technique has recently been used to measure the proper motion of a lens 
seen toward the Small Magellanic Cloud (SMC) and so demonstrate that the lens 
almost certainly resides in the SMC (Afonso et al. 1998; Albrow et al. 1999; Alcock 
et al. 1999; Udlaski et al. 1998; Rhie et al. 1999). However, such caustic crossing 
events are rare and the great majority of them are binaries (and hence may not be 
representative of the lens population as a whole). 

A second type of information can come from parallax measurements. If the 
event is sufficiently long, then the normal light curve is distorted by the accelerated 
motion of the Earth about the Sun allowing one to measure the Einstein radius 
projected onto the observer plane, r E , 

r E = D9 E , (3) 

Gould (1992). Several parallaxes have been measured for bulge events (Alcock 
et al. 1995; private communication D. Bennett 1998), but all for events that 
are substantially longer than typical. It would be possible to routinely measure 
the parallaxes of microlensing events by launching a satellite into solar orbit 
(Refsdal 1966; Gould 1995b). The event would have a different time of maximum 



magnification, to, and different impact parameter, /3, as seen from the Earth and the 
satellite. From the differences in these quantities, At and A/3, and using the known 
Earth-satellite separation, <i©_ s , and known angle 7 between the line of sight and the 
Earth-satellite vector, one could reconstruct both the size of the projected Einstein 
ring, te, and the direction of motion, 0, relative to the satellite-Earth vector, 

- d e - s \ sin -f\ a_ A P u\ 

rE -[(Apy + (At /t E y]^ tan0 "A^ (4) 

A rather technical but in the present context very important point is that it 
is significantly easier to measure the difference in times of maximum, Ato, than it 
is to measure the difference in impact parameters, A/3. There are two interrelated 
reasons for this which are investigated in detail by Gould (1994b; 1995b), Boutreux 
& Gould (1996), and Gaudi & Gould (1997). First, the sign of the impact parameter, 
/3, measured by a single observer is intrinsically ambiguous because the light curve 
contains no information about the side of the lens on which the source passes. 
Hence, from the two individual impact-parameter measurements, /3® and /3 S , it is 
possible to reconstruct four different values of A/3 = ±(/3 e ± f3 s ). Second, one must 
determine (3 and t from the light curve simultaneously with three other parameters, 
tE, Fq, and Fb, the latter two being the fluxes from the source star and any unlensed 
background light that is blended into the photometric aperture of the source. While 
to is virtually uncorrelated with any of these three other parameters, /3 is highly 
correlated with all of them, in particular with Fb. Hence (3 (and so A/3) is more 
poorly measured than to (and Ato). While it is possible to break the four- fold 
discrete degeneracy, this requires measurement of a higher-order effect. 

No dedicated parallax satellite is currently planned. However, the Space 
Infrared Telescope Facility (SIRTF) could be used to measure parallaxes of at 
least some events. Because SIRTF makes its measurements in passbands that are 
inaccessible from the ground, the relative blending between the Earth and satellite is 
completely unconstrained, so measurement of A/3 is not simply difficult, it is virtually 
impossible. Nevertheless, if Ato is well constrained by Earth-satellite observations, 
then it is possible to determine A/3 from vigorous ground-based observations (Gould 
1999a). 

In the best of all possible worlds, one would measure both #e and te- These 
(together with the routinely measured tE and the approximately known source 
distance) would then yield a complete solution for M, _D ol , and v. (e.g. Gould 



1995c). For example, 



c 2 



M = — e E r E . (5) 

At present, this is possible by ground-based measurements only for certain rare 
classes of events (Hardy & Walker 1995; Gould & Andronov 1999; Gould 1997). 
If there were a parallax satellite, then it would also be possible for those rare 
events which happened to be accessible to proper- motion measurement. However, 
astrometric microlensing opens the possibility, at least in principle, that such 
complete measurements might be made for a large unbiased sample of events in the 
future. 



2. Astrometric Microlensing: Promise and Limitations 

As Boden et al. (1998) discuss, astrometric measurements are sensitive to two 
distinct effects. First, the center of lensed light from the source is displaced from the 
actual position of the source by, 

where u = fi(t — £o)/#e is the projected position of the source relative to the lens 
in units of angular Einstein radius assuming rectilinear motion as would be observed 
from the Sun. That is, 
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(7) 



This deviation traces out an ellipse with semi-major and semi-minor axes, 

1 (3 

9a = w+W 9ei 9b = W+z) 9e ' (8) 

The major axis is aligned with the direction of motion of the lens relative to the 
source. Hence, by measuring this effect, one can solve for both 0-& and the direction 
of motion. A measurement of #e is often called a "proper motion" measurement 
because, from equation (|1|) it can be combined with the known Einstein time scale 
to yield the magnitude of the proper motion. However, in the case of astrometric 
measurements, it also yields the direction, 0, and so the full vector proper motion, fi. 

The second effect is a parallax deviation caused by motion of the Earth about 
Sun. The exact formula for the combined parallax and proper motion effect can be 



found by substituting 

n x n x a 

^ = u (9) 

r E 

into equation (g). Here n is the unit vector in the direction of the source, and a is 

the position vector of the Earth relative to the Sun. Thus, the magnitude of the 

perturbative term is ~ AU/fs, which might be ~ 10%-30% for typical lensing events. 

This would seem to imply that one could determine the parallax (fg) about 10% to 

30% as accurately as the proper motion (#e) for the same set of measurements. 

Unfortunately, the situation is not so favorable. The perturbation in equation 
([J) is not directly observable because there are no comparison observations from 
the Sun. Consider the limit £e "C Jt/2tt which is typical for events seen toward the 
Galactic bulge. In this case, the Earth's velocity would barely change during the 
event or even for the first few £e after it. One would then see the same ellipse as 
described by equation (H), but with a different (3 from the one that would have been 
seen from the Sun. Ninety days after the event, the direction of apparent source 
motion would have changed by an angle v^/v where v = te/^e is the speed of the 
lens projected onto the observer plane, and % = 30kms _1 is the speed of the Earth. 
According to equation (|^), this would introduce an astrometric deviation of order 
(f®/v)(tE/90day)6 , E. Since 9e, is essentially determined from the major axis of the 
ellipse (eq. [5]), which is approximately 2 _1 / 2 6*e, the relative size of the astrometric 
parallax and proper-motion measurement errors is roughly given by 



_ Q'fE.ast/^E v 60 da y 

ast/ast — I n r ^ J j. 

&6 E ,ast/VE "% £e 



1 ast/ast = 77T ~ 7 • (,-LUJ 



In fact, equation (10) is too optimistic in that it implicitly assumes that the direction 
of the major axis of the ellipse can be determined with infinite precision. As we will 
show in § 5, this is very far from the case. Hence, the true ratio of errors is generally 



larger than implied by equation (10). For typical bulge-bulge lensing events, 

v ~ 800kms _1 . For typical halo events seen toward the LMC, v ~ 300kms _1 . For 

lenses in the LMC, v is a factor 3 to 10 higher yet. 

The above analysis implies that astrometric microlens parallax measurements 
are several orders of magnitude less accurate than proper-motion measurements. This 
would not present much of a problem if very accurate proper motion measurements 
could be made with a modest amount of observing time. However, as we will show 
in § 5, even for bright (I ~ 15) sources seen toward the bulge, proper-motion 
measurements accurate to 5% require about 5 hours of observations. Very few events 



seen toward the LMC are brighter than V ~ 20 and therefore ~ 40 times more 
observing time is needed to achieve the same precision. Hence, this analysis appears 
to imply that no more than a few accurate microlens parallaxes could be obtained 
in any reasonable observing program. 



3. Photometric microlensing parallaxes with SIM 

SIM is not designed to do photometry, and it would seem completely hairbrained 
to waste this precision astrometric instrument on measurements that could be done 
more efficiently from the ground using telescopes with collecting areas that are 
several orders of magnitude larger. Nevertheless, two unrelated factors combine to 
make SIM the ideal device to measure microlens parallaxes photometrically (rather 
than astrometrically) . 



3.1. Photometry with SIM 

First, SIM works by counting photons as a function of position in the interference 
pattern in order to find the centroid of the central fringe. The photons are distributed 
in this fringe as NF(8)d6, where N is the total number of photons in the central 
fringe, 

^>=4 cos2 (4> "'-is- < n > 

d is the distance between the mirrors, and A is the wavelength of the light. The 
astrometric precision is given by, (e.g. Gould 1995a), 

r /■ /<71n F\ 2 i - 1 / 2 

a e = N- 1 / 2 [J d9F(9) (^f) ] = N~V%, (12) 

and hence the fractional photometric precision (u p h = N~ x ' 2 ) is related to the 
astrometric precision by 

^ P h = /■ (13) 



3.2. One-dimensional photometric parallaxes 

Second, as discussed in the introduction, the real problem in obtaining 
microlens parallaxes photometrically is that the microlens parallaxes are inherently 
two-dimensional. In effect, by measuring Aio, one determines cos0/fE, and by 
measuring A/3, one determines sin (f)/r E , where is the angle of source-lens relative 
motion with respect to the direction of SIM-Earth axis at the moment when the 
event is a maximum as seen from the SIM-Earth midpoint. It is only by measuring 
both of these quantities that one can determine r E . Since A/3 is difficult to measure, 
obtaining a precise r E is also difficult. 

As discussed in the introduction, it is possible in principle to break the 
degeneracy in Af3 photometrically if the photometry is good enough. We will show 
in § 5 that SIM photometry is sufficiently precise for this task provided that the 
observations are carefully planned. However, it is also the case that SIM astrometric 
measurements by themselves often determine (from the orientation of the ellipse) 
with sufficient precision to break the degeneracy in Af3. In these cases, r E can be 
determined from a measurement of At alone. In this paper, we will consider both 
methods of breaking the degeneracy, but in the remainder of this section we will 
assume that the degeneracy is broken astrometrically. This will allow us to estimate 
r p h/ast> the relative precision of SIM photometric parallax measurements to SIM 
astrometric proper-motion measurements. 

As currently designed, SIM will fly in a SIRTF-like orbit, drifting away from the 
Earth at about 0.1 AU per year. Let the distance at the time of the observations be 
d®_ s . Then r E can be determined from the measured At (and the known value of 

4>) by, 

t E 
r E = d e _ s — — |cos0|. (14) 

Gould (1999a) analyzed how to optimize measurements of Ato when he investigated 
microlens parallaxes with SIRTF. For photon-limited photometry, one should 
concentrate the measurements near times t± before and after the peak, where 
t± = to ± (5/3) 1 / 2 / 9tE- The error in to,s is then approximatley given by 

<7«,,. ~ ™ -1/ V ph ff* (/?£0.5), (15) 

where n is the total number of measurements in these two regions. We will assume 
that the ground-based measurements to determine to,® have similar precision as the 



SIM measurements. Equation ([14]) then implies that the fractional error in te is 
given by, 

^E = 21/2^, = 2 1/2^| gec0 |_ (16) 

te |At | a©- s 

Then, assuming that a similar number of measurements are used to determine the 
semi-major axis of the ellipse, 9 a , and to determine At (they are somewhat the 
same measurements), the ratio of the photometric precision of the parallax to the 
astrometric precision of the proper motion is 



CTfy/rE r E a p h(/3/O.3)|sec0 |/(j @ _ 



J- ph/ast = — — Yq- ~ 777i /oN , (.1 < J 



where we have approximated equation (|j) as # a = 6*e/3. Using equations @ and 
( PI) , this can be rewritten 

rph / ast ~^^^ |seC0l ^^O^Melo^AuJ |S6C01 ' (18) 

where we have adopted 9f = 2.5 mas, which is appropriate if the flux-weighted 
harmonic mean wavelength of the source is 0.8 /mi, and the mirrors are separated by 
10 m. Equation fll8f ) implies that for typical lenses, the photometric parallax will be 
of comparable precision to the astrometric proper motion in sharp constrast to the 
large ratio for astrometric parallaxes found in equation (|i~0|) . 

In fact, the | sec0| dependence in equation ( p|) is too pessimistic because we 
have ignored all photometric information about A/3. We show in § 5 that except 
for the case cos0 ~ (where the discrete degeneracy becomes astrometrically 
incorrigible) it is possible to essentially eliminate the | sec0| term in equation ( |T7| ) 
using a combination of astrometric and photometric data. 



4. Simulations 

The estimates given in the previous two sections are useful because they 
elucidate the relation between the physics of the event and the measurement process 
on the one hand and the errors in the microlensing parameters on the other. By the 
same token, however, they cannot capture the full range of experimental conditions, 
and so are necessarily approximate. The actual errors for any given event will depend 
both on the precise event parameters and on the observational strategy. While a 
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full investigation of the best observational strategy lies well beyond the scope of the 
present study it is important to make a rigorous calculation of the statistical errors 
for some representative examples in order to obtain more precise estimates and to 
investigate more subtle effects that are not captured by the rough analysis given 
above. 

For this purpose, we consider a set of somewhat idealized observations. First, 
we assume that the principal measurements are carried out at uniform time intervals 
that are short compared to £e, beginning when the magnification first reaches 
A = 1.5 and ending at a time that will be determined below from signal-to-noise 
(S/N) considerations. This is quite reasonable for observations toward the LMC 
(near the ecliptic pole) but is obviously not really possible toward the bulge (near the 
ecliptic). For the bulge, we therefore assume that the measurements are interrupted 
when the bulge is within 60° of the Sun. This measurement strategy can actually be 
very far from ideal, and we modify it somewhat in § 4.2, below. Second, we assume 
that the exposure times are of equal durations, so that the S/N is better near the 
peak of the event. Third, we assume 9f = 2.5 mas and rf ffi _ s = 0.2 AU, although the 
first will clearly vary from star to star, and the second will change during the course 
of the mission. For the LMC (at the ecliptic pole), the time of year at which the 
event is discovered plays no role, but for the bulge it does. We will assume that the 
bulge field lies at B = — 6° from the ecliptic, close to the (northern) winter solstice. 
The Earth-satellite separation projected onto the sky is therefore at a maximum 
at the summer solstice and is d e _ s |(l — cos 2 ip cos 2 B)] 1 I 2 ~ d®- s \ sin ip\ at other 
times of the year, where ip is the phase of the Earth's orbit relative to the autumnal 
equinox. As discussed in detail by Gaudi & Gould (1997), the orbital phase if) 
has two conflicting effects. First, when cosip ~ (near the summer solstice), the 
SIM-Earth projected separation is a maximum, and hence the measurement errors 
of Ato and A/3 are reduced to a minimum. On the other hand, the relative velocity 
of SIM and the Earth projected onto the plane of the sky is at a minimum, and 
breaking the degeneracy in A/3 depends critically on this relative velocity. Hence, it 
is most difficult to break the degeneracy at the summer solstice. As the phase moves 
away from the summer solstice, the projected separation slowly declines (making 
the measurements of Ato and A/3 less accurate) but the projected velocity difference 
rapidly increases (allowing more secure degeneracy breaking). If one relies on 
photometry to break the degeneracy, the optimal events are those that peak about 
45 days from the summer solstice (see e.g. Fig. 6 from Gaudi & Gould 1997). On 
the other hand, if it is possible to break the degeneracy astrometrically, then events 
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that peak at the solstice are optimal since the measurement errors are reduced by 
~ 2 -1 / 2 . In this paper, we will primarily simulate events peaking at ip = 225°, i.e. 
about May 7. However, we will also discuss events that peak at the summer solstice 

(V> = 270°). 

Next, we will assume for definiteness that the ground-based photometric 
observations have the same precision as the SIM observations. Finally, we will 
ignore blending in the SIM measurements, except blending by the lens. Blending 
will have an important impact on the overall precision and hence the strategy for 
SIM measurements but as we show below, it will not affect the main conclusions of 
this paper which concern the relative precision of SIM astrometric and photometric 
microlens parallaxes. This is especially so toward the bulge which will be the main 
focus of analysis. A proper treatment of blending would therefore make the paper 
substantially more complex without clarifying any of the central points. Hence, we 
defer consideration of this important effect to a later paper on observational strategy. 

Why can blending be ignored to first order in this analysis? First, Han & 
Kim (1999) have shown that all potential blends lying more than 10 mas from the 
source can be eliminated by the SIM observations themselves. Since the density 
of field stars having even a modest fraction of the source flux is much less than 
10 4 arcsec -2 , this essentially eliminates all blends not directly associated with the 
event, namely the lens itself, binary companions to the lens, and binary companions 
to the source. For bulge events, 10 mas corresponds to 80 AU, so a substantial 
fraction of binary companions are also eliminated. Second, to minimize observation 
time, SIM observations must be almost entirely restricted to very bright stars. This 
means clump giants toward the bulge and either clump giants or early main-sequence 
stars toward the LMC. The chance that a bulge clump star has a companion within 
80 AU and with more than a few percent of its own flux is small because their 
progenitors are about 50 times fainter than the stars themselves. The primary 
effect of a few percent blend would be to change the shape and orientation of the 
proper-motion ellipse and (assuming the shape change went undetected) to therefore 
change the inferred direction of the lens-source relative proper motion, also by a few 
percent. This would in turn affect the parallax inferred from At which depends 
on this direction through the angle 0. See equation (|l4j) . However, this effect 
will also be a few percent. As we will show in § 5, it is quite possible to achieve 
accuracies of a few percent for bulge events, so a careful investigation of the effect 
of blending on parallax and proper-motion measurements should be undertaken as 
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part of a more thorough analysis of the problem. Unfortunately, a proper analysis of 
blending from binary companions to the source requires simulated fits to the entire 
diffraction pattern, not just the centroid, and so is substantially more involved than 
the present study. By contrast, low- level blending by the lens can be treated within 
the framework of the centroid analysis given here and we therefore include it. 

The situation is more complicated toward the LMC because the chance that 
an early main sequence star has a companion of comparable brightness is larger, 
probably a few tens of percent. Even clump stars have brighter companions toward 
the LMC than toward the bulge because they are younger and so have brighter 
progenitors. Also, the 10 mas limit on detecting blends directly (Han & Kim 
1999) corresponds to 500 AU toward the LMC compared to 80 AU toward the 
bulge. Nevertheless, even toward the LMC, the majority of sources will not have 
companions with more than 10% of the source flux and hence even here it is 
appropriate to ignore blending by companions in a first treatment. 

Toward both the bulge and the LMC, it is very unlikely that the lens itself will 
contribute more than a small fraction of the source light if the source is bright. We 
therefore allow for lens blending in our simulated fits but assume that the actual 
blending is very small. 

Note that we will not ignore blending in the ground-based photometric 
observations, since there is no way to eliminate field star blends for the ground-based 
observations. 



4.1. Parameterization 

We will simulate simultaneous observations from SIM and from the ground. 
There will be four measured quantities: 1) G 1 , the flux observed from the ground, 
2) G 2 , the flux observed from SIM, 3) G 3 , the x astrometric position, and 4) G 4 , the 
y astrometric position. These give rise to four observational equations, 

G\t) = F}A( U(B ) + F b \ G 2 (t) = F 2 A{u SIM ) + F 2 (19) 

and 

[G 3 (t), G\t)} = " SI ^ 6 E - asiuKs + H,t + o -9 b (20) 
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where A(u) = (u 2 + 2)/[u(u 2 + 4) 1 / 2 ] is the magnification, /j, s is the absolute proper 
motion of the source, 0q is the true position of the source at time t = 0, 

n x n x a . . 

° = — Sij-- (22) 

n s is the parallax of the source, and Usim and u e are defined similarly to equation 

(|), i.e. 

AU 
usim = u + Ka S iM, u e = u Q + «a®, k = -^, (23) 

with 

u = (r cos — (3 Q sin 0, r Q sin + /3 cos 0), r Q = — — — . (24) 

*E,0 

The terms in equation ( p0|) can be understood as follows. The first term is the 
ellipse characterized by equation (|8]), modified by the motion of the Earth (i.e. 
u — > usim)- The second term is the parallactic motion of the source. The third 
and fourth terms represent the ordinary proper motion and position of the source. 
The last term is the perturbation due to the luminosity of the lens which is written 
out explicitly in ([H]). Here the first term is the difference between the "ellipse" in 
equation ( p0|) and the relative lens-source position (— Usim9&), while the second is 
the relative parallax of the lens and the source. 

There are a total of 15 parameters: £ ,©) Ad, and £e,© are the standard event 
parameters as it would be seen from the Sun, is the direction of source motion 
relative to the lens with respect to the SIM-Earth direction, Be is the Einstein radius, 
K is the inverse projected Einstein radius (normalized in AU), ir s is the parallax of 
the source, /j, s is the proper motion of the source, 9 is its position at t = 0, F} and 
F 2 are the source fluxes as received by the Earth and satellite observatories, and F^ 
and F 2 are the background fluxes. 

To determine the uncertainties in these parameters, we evaluate the covariance 
matrix c^ (e.g. Gould & Welch 1996) 

L-l L Y^Y^ -2^ (tk) 9G (t k ) , s 

c = b • 6 « = §1> H a* da, ■ (25) 

where ai ... a±5 are the fifteen parameters, tk are the times of the observations, and 
G]~i is the error in the measurement of G l at time tk- We enforce the condition of 
weak blending by setting F^ = F 2 = after taking the derivatives in equation (p5|). 
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4.2. Modification of observing strategy 

The full problem of optimization of SIM microlensing observations lies outside 
the scope of this paper, but it is straight forward to determine the optimum duration 
of observations once the (arbitrary) strategy of uniform observations has been 
adopted: one changes the interval over which the observations are carried while 
holding the total observing time fixed and inspects the resulting errors. We carry out 
this exercise and find that the optimal duration to determine r& is short, typically a 
few tens of days for various combinations of parameters, while the optimal duration 
to determine #e is well over 100 days. The reason for this is clear. The measurement 
of Te is determined primarily from photometry, and the photometric microlensing 
event is essentially over after 2£e- On the other hand #e is determined from the 
astrometric event which lasts many £e- Only after the astro metric event is essentially 
over is it possible to determine /j, s and so remove the correlation between this 
parameter and 9e- 

We address this inconsistency of time scales by modifying the observational 
strategy. We take observations uniformly over various intervals but with fixed total 
observing time, T, and then take three additional observations, each with observing 
time T/20, at 3 months, 9 months, and 12 months after the peak of the event. Thus, 
the total observing time is 1.15 T. We then find that the precisions of both f-& and 
#e are roughly constant when the continuous observations last anywhere from 30 to 
120 days. For definiteness, we adopt 50 days uniformly. 



4.3. S/N: assumptions and scalings 

We assume that all sources are / = 15 which is typical of the brighter 
microlensing events seen toward the bulge but is much brighter than any event seen 
so far toward the LMC. In our calculations, we normalize the astrometric precision 
by assuming that 4/zas accuracy (in one dimension) can be achieved in 1 minute of 
observation on an I = 11 star. That is, our fiducial / = 15 stars require 40 minutes 
to reach 4/ias. We allow a total of 5 hours of observation for each event. It is then 
straight forward to scale our results to other assumed conditions. For example, 
for an / = 20 source, the errors reported in the next section must be multiplied 
by 10. Alternatively, the same errors could be achieved by allowing 500 hours of 
observations. If our astrometric error estimate proves too pessimistic, so that it is 
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possible to achieve 4/xas precision in a minute on an / = 12 star, then the errors 
should be divided by 1.6. 



5. Results 

5.1. Bulge 

We consider two geometries. The first is a bulge line of sight at 6° from 
the ecliptic. The source and lens are both in the bulge, with D os = 8kpc and 
D Q \ = 6kpc. Hence D = 24kpc. (See eq. ||.) The speed of the lens relative to the 
observer-source line of sight is v = 200 km s" 1 . We vary M, <b, and (3. Formally 
(3 is defined as the impact parameter of the event as seen from an observer at the 
Earth-SIM midpoint but in practice it is very similar to the /3 observed from the 
Earth. All events are assumed to peak (as seen from the midpoint) on May 7, i.e., 
45 days before the (northern) summer solstice. This is the most favorable time to 
break the discrete degeneracy in A/3 photometrically (see Gaudi & Gould 1997), but 
the intrinsic errors in A/3 and Ato are larger by ~ 2 1 / 2 than they would be at the 
summer solstice. We will therefore later investigate whether the discrete degeneracy 
in A/3 can be broken astrometrically so that observations could take place at the 
solstice. In order to better understand this and several other issues, we conduct 
two sets of simulations, one using SIM measurements alone (both astrometric and 
photometric) and the other combining astrometry with photometry from both SIM 
and the ground. 

Table 1 shows the results. The first three columns are the input parameters 
M, f3, and <fi. Columns 4 and 5 show the errors in cb from SIM measurements only 
and with the addition of ground-based measurements. Columns 6 through 11 show 
the fractional errors in various quantities both without and with ground-based 
photometry. Columns 6 and 7 show the errors in fg, columns 8 and 9 show the 
errors in 9e, and columns 10 and 11 the errors in M = f-^O-^c 2 / 'AG . Finally, columns 
12 and 13 show the fractional errors in D = te/Oe, and in n s based on combined 
data from SIM and the ground. 

These results are in rough qualtitative agreement with the predictions of 
equations (|T0D and ([181) . In particular, they confirm that the fractional error in fg 
is orders of magnitude larger than the fractional error in #e if one is restricted to 
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SIM data, while the two errors are comparable if one combines astrometric with 
photometric data. 

However, there are a number of additional important conclusions that can be 
drawn from Table 1. First, restricting consideration to impact parameters (3 < 0.5, 
the fractional errors in M, D, and ir s are all typically about 5% although these errors 
do vary somewhat in particular cases. This means 5 hours of observation produce 
very precise individual solutions for the mass, distance, and velocity of the lens, and 
also for the distance and velocity of the source, implying that a few hundred hours 
of SIM time could yield a very detailed inventory of the material between the Sun 
and the Galactic center. 

Second, while the errors in r~E do deteriorate toward = 0, the trend is not 
as drastic as predicted by equation ([T8fl. The fundamental reason for this is that 
the continuous degeneracy in A/3 is not very severe, so that if one assumes that the 
discrete degeneracy is broken, then there is actually quite a lot of information about 
this component of r^ in the photometric measurements. In fact, the values in Table 
1 implicitly do assume that the discrete degeneracy is broken. This is because they 
are based on equation ( p5| ) which gives a purely local error analysis. 

Recall that there are two discrete degeneracies. We focus initially on the one 
that affects the magnitude of A/? and defer consideration of the one that affects 
only the sign of A/3. To determine whether this discrete degeneracy is broken 
photometrically, we examine the work of Gaudi & Gould (1997), in particular 
their Figure 6. Under the observational conditions they considered, the discrete 
degeneracy is broken 90% of the time for M = 0.3 M , tp = 225°, and satellite 
separations d®- s = 0.4 AU. This is twice the separation that we have assumed. 
However, Gaudi & Gould (1997) have assumed photometry errors of 1% for the 
Earth and 2% for the satellite for a total of about 70 observations. If our 5 hours 
of observing time were divided among 70 observations, the photometric precision 
would be 1.2%. That is, our assumed effective errors are smaller by a factor 
[(1.2 2 + 1.2 2 )/(1 2 + 2 2 )] 1/2 = 0.75. For small Earth-satellite separations, there is a 
direct trade off between measurement error and satellite separation, so our 0.2 AU 
separation corresponds to 0.25 AU in their simulations. Inspection of the Gaudi 
& Gould (1997) Figure 6 shows that the degeneracies would be broken 70% of the 
time. To determine how this effectiveness scales with lens mass, we turn to Figure 
4 from Gaudi & Gould (1997). This shows that degeneracy breaking becomes more 
difficult at lower masses. The figure is drawn for the case <i©_ s = 1 AU whereas the 
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argument just given implies that with 5 hours of observation, our 0.2 AU separation 
is equivalent to 0.25 AU in the Gaudi & Gould (1997) simulations. Hence, comparing 
Gaudi & Gould (1997) Figures 4 and 6, we estimate that 35% of the degeneracies 
would be broken for M — 0.1 M & . If the exposure times were multiplied by a factor 
of 4 to 20 hours, this fraction would rise to about 70%. 

We now turn to the question of how well the degeneracies can be broken 
astrometrically. For the geometry considered here, f-& = 7.5(M/O.3M ) AU. Hence 
A/3 ~ rf ffi _ s /fE ~ 0.025, which is quite small compared to typical values of (3. This 
means that for most events the discrete degeneracy will be between solutions with 
A(3 ~ 0.025 and A(3 ~ 2(3 ~ 0.5. Since tan0 = A/3/Ai , the high A(3 solution 
will almost always correspond to angles 90° ± 2° or 270° ± 2°, while the low A/3 
solution (assuming that it is the real one) will be at some random angle. Thus to 
distinguish the two solutions, one must have independent information about <fi with 
errors that are a factor ~ 3 smaller than |90° — <p\ or |270° — <p\. Is this well 
constrained by the observations? Looking at column 4 of Table 1, seems to be 
very well constrained. However, this precision measurement is based primarily on 
the photometric measurements of A/3 and Ato and so implicitly assumes that the 
discrete degeneracy has been broken. Hence we should use only SIM data (column 
5). 

We see from column 5 of Table 1 that the errors in <fi are generally small for 
(f) < 60° but deteriorate toward <fi = 90°. This means that the discrete degeneracy is 
broken astrometrically at the 3 a level for <p < 60° but cannot be broken if the angle 
gets close enough to 90°. Table 1 does not have sufficient resolution to determine 
the transition but we find by more detailed calculations that for (3 = (0.2, 0.4) this 
occurs at (70°,65°), for M = 0.1 M , at (75°,70°), for M = 0.3 M , and at (80°, 75°), 
for M = 0.5 M . Hence, the degeneracy is usually broken astrometrically, but less 
frequently at low masses. Since the degeneracy is more difficult to break at low 
masses both photometrically and astrometrically, it would be prudent to commit 
more observation time (say 20 hours rather than 4 hours) to the shortest events 
(which are likely to be low mass). 

The fact that the degeneracy can be broken astrometrically for most events 
seems to argue against restricting observations to periods that are 45 days from 
the summer solstice. Recall that we adopted this restriction in order to permit 
better photometric degeneracy breaking which now no longer seems so necessary. 
However, we find that for events peaking at the solstice, the errors in (when the 



Earth-based observations are ignored) are substantially higher than the values in 
column 5, implying that it is not often not possible to break the discrete degeneracy 
astrometrically or photometrically at the solstice. 

As we have noted, the above discussion actually applies to only one of two 
discrete degeneracies, the one involving two different magnitudes of Aj3. This is the 
more important degeneracy because it affects the estimate of the size of fg and so of 
the mass, distance, and speed of the lens. However, there is also another degeneracy 
involving the sign of A/3 but not its magnitude. Boutreux & Gould (1996) and 
Gaudi & Gould (1997) refer to the first of these as the "speed degeneracy". We call 
the second the "direction degeneracy". The direction degeneracy becomes difficult 
to break when \Aj3\ <C |A£o|/£e, i-e., when ~ 0° or ~ 180°, so that the two 
degenerate solutions are close in 0. From column 5 in Table 1, we find that the error 
in from the astrometric data alone is generally quite small for = 0, and hence 
is adequate to break the direction degeneracy unless lies within a few degrees of 
either 0° or 180°. However, in this case the effect of the degeneracy is very small. 

Finally, we note that we have examined the errors in F^/ F s (although we do 
not display them). We find that they are typically a few percent, implying that a 
luminous lens could be detected if it were more than a few percent of the flux of the 
source. 



5.2. LMC 

The second line of sight is toward the LMC at the south ecliptic pole. The source 
lies at D os = 50kpc, while the lens is assumed to be in the halo at D \ = 15kpc. 
Hence D = 21 kpc, which is very similar to the bulge value. This means that at fixed 
mass, the bulge events considered in § 5.1 will have about the same te as the LMC 
events considered here. The speed of the lens relative to the observer-source line of 
sight is v — 250 km s -1 , slightly larger than for the bulge. We have again assumed 
an I — 15 source to allow easy comparison with results from the bulge. We discuss 
the implications of this highly unrealistic assumption below. 

Table 2 shows the results. They are qualitatively similar to those for the bulge. 
The largest difference is that the fractional error in ir s is larger which simply reflects 
the fact that the LMC is more distant. 
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Of the eight microlensing events detected by Alcock et al. (1997b) during their 
first two years of observations, none of the sources were brighter than V = 20 (after 
removing blended light) which is the nominal limit for SIM. Future microlensing 
surveys could improve the rate of detection by an order of magnitude (Gould 1999b; 
Stubbs 1998). However, only a factor of three of this improvement would be due to 
the coverage of a larger area. The rest would come from going deeper which would 
not yield any more bright sources. Hence, the total rate of events that are accessible 
to SIM will not be high. Most of the usuable events that are detected are likely to be 
close to the magnitude limit. For astro metric measurements V = 20 is equivalent to 
/ = 19 because at fixed magnitude there are fewer photons in I and they are longer. 
Hence, for the same total observing time (5 hours), the errors will be about a factor 
of 40 1 / 2 larger than listed in Table 2. That is, the mass and distance estimates will 
be accurate to about 20 to 40%. This would be an acceptable level of precision to 
resolve the question of the nature of the lenses assuming that more than a handful of 
events can be measured. The errors in the measurement of <fi without making use of 
the Earth-based observations (column 5) are typically 15° (after including the factor 
40 1 ' 2 ). Hence in many cases it will not be possible to break the A/5 degeneracy 
astrometrically. 

To determine whether the photometry is sufficiently precise to break the 
degeneracy, we compare our simulation with that of Boutreux & Gould (1996) who 
specifically considered an Earth-satellite separation of 0.26 AU which is close to our 
value of d e „ s = 0.2 AU. In their Monte Carlo simulation, the "speed degeneracy" 
(between different scalar values of A(3) was broken 40-60% of the time in the mass 
range 0.1-1 M Q . We find that their assumed photometric errors are about equal to 
the ones assumed here (for a total of 5 hours of observation at V — 20). Therefore, 
it seems likely that the SIM photometric observations would be adequate to break 
this degeneracy in the majority of cases. 

Inspection of column 5 from Table 2 shows that the direction degeneracy will 
usually be broken astrometrically. The simulations of Boutreux & Gould (1996) 
show that it is about equally difficult to break the direction degeneracy as the 
speed degeneracy. Thus, it should usually also be possible to break this degeneracy 
photometrically. In any event, as in the case of the bulge, the direction degeneracy 
is much less important than the speed degeneracy. 

While 5 hours of observation is adequate to characterize events toward the LMC, 
it would seem advisable to invest more time for these events. First, there will be very 
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few of them so the cost in total observing time is not high. Second, if the events are 
in the halo, one would like to know their mass more accurately than 30% if possible. 
Finally, if the lenses reside in the LMC, then 5 hour exposures would yield only a 
non-detection of parallax rather than a measurement. While this (together with the 
proper-motion measurement) would be highly significant in demonstrating that the 
lenses were not in the halo, it would provide only limited clues as to their nature. In 
brief, the scientific return of extra observation time would be high for LMC events. 
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